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Inspired by the newly observed Dsj (3040)^ state, in this work we systemically study the two- 
body strong decays of P-wave charmed-strange mesons with the first radial excitation. Under the 
assignment of = I"*"), i.e. the first radial excitation of _Dsi(2460)^, we find that the width 

of _Dsj(3040)^ is close to the lower limit of the BaBar measurement. This indicates that it is 
reasonable to interpret _Dsj(3040)^ as the first radial excitation of -Dsi(2460)^. Our calculation 
further predicts that 0~ + 1~ channels e.g. D'^ K*'^ , D'^K*'^ and D^cfi are important for the search 
for _Dsj(3040)"'". To help future experiments finding the remaining three P-wave charmed-strange 
mesons with the first radial excitation, we present the predictions for the strong decays of these 
three P-wave charmed-strange mesons. 

PACS numbers: 13.25.Ft, 12.39-x 



I. INTRODUCTION 

With the new observation of Dsj meson, the spectrum of the charmed-strange state is becoming abundant. So 
far, there exist six established charmed-strange mesons £'s(1968), i:'*(2112), D*q{2317), 1551(2460), £'si(2536), 
D*2(2573) listed in Particle Data Group (PDG) Q, which can be categorized as three doublets in terms of 
the heavy quark limit: H (O^,:-) = (L»,.(1968), D*(2112)), S = (0+,l+) = (i:i:o(2317), i:»3i(2460)) and 
T = (l+,2+) = (D,i (2536), £1*2(2573)). Two years ago, a new charmed-strange meson D;i(2710) with = 1" 
was firstly announced by the Babar Collaboration Q and confirmed by the Belle Collaboration later Very re- 
cently the Babar experiment found D*^{2710) again in the D*K invariant mass spectrum Q. Another newly observed 
charmed-strange meson is D*j(2860), which was observed in both DK j^land D*K channels JSl. The phenomeno- 
logical proposals of the quantum number of £'*j(2860) include = 3~ (1, [B| and = 0+ [aJ0, Hi- As indicated 
by the Babar experiment, = 0^ assignment for D*j(2860) is forbidden according to the parity conservation since 
the D*K decay mode of L>*j(2860) was observed in Ref. Q. A series of theoretical work 0, [1, 0, S, i, E HH, E [11 
relevant to i:>*i(2710) and i:»;j(2860) were carried out. 

Besides the observations of D*i(2710) and D*j(2860) by analyzing the D*K invariant mass spectrum in inclusive 
e+e" interactions 3, Babar also announced a new charmed-strange state £'s,/(3040) with the mass M = 3044 ± 
8(stat)tf (syst) MeV and the width T = 239 ± 35(stat)t4^(syst) MeV Q. The observation of £)sj(3040) not only 
makes the spectrum of the charmed-strange meson abundant (the mass spectrum of the observed charmed-strange 
mesons are listed in Fig. [T|), but also stimulates our interest in exploring its underlying structure. 

As indicated by the Babar Collaboration, Z?s,/(3040) was only observed in D*K channel while not find in DK 
decay mode. Thus, its possible quantum number includes = 1+, 0~, 2~, • • • . Since Dsi{2710){J^ — 1~) is the 
first radial excitation of L'*(2112) and the mass of Dsj{i04:0) is far larger than that of Z?si(2710), thus we further 
exclude 0~ assignment, viz. the first radial excitation of Ds{1968) for Dsj{304:0). In Ref. Matsuki, Morii and 

Sudoh once predicted the mass of cs state with ■nP'^^^Lj = 2^ Pi. m = 3082 MeV, which is close to the experimental 
value of the mass of Z?sj(3040). Thus, 1+ assignment to Z)sj(3040), the first radial excitation of Dsi{2AQQ), becomes 
the most possible. 

If Dsj(3040) as the radial excitation of P-wave charmed-strange state is true, further experiment is of the potential 
to search the rest three radial excitations of P-wave charmed-strange states. Thus, a systematical phenomenological 
study of the strong decay mode of P-wave charmed-strange mesons with the first radial excitation is an important and 
interesting topic. By this study, we will not only obtain the information the decays of these P-wave charmed-strange 
mesons, but also can test 1+ quantum number assignment to Z?sj(3040) comparing calculated decay width with the 
experimental data. 
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In this work, we will be dedicated to the study of the strong decay modes of P-wave charmed-strange mesons with 
the radial excitation by the ^Pq model [H, [H, O [H, [H, HO, EH . Further we will obtain the information of the order 
of magnitude of the strong decay modes of Dsj(3040). 
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FIG. 1: The mass spectrum of the observed charmed-strange mesons and the corresponding strong decay modes observed in 
experiment 



The paper is organized as follows. After the introduction, we briefly review the ^Pq model and present the formula- 
tion of the strong decays of P-wave charmed-strange mesons with the radial excitation. Finally, the numerical result 
will be shown. The last section is a short summary. 



II. THE STRONG DECAY OF P-WAVE CHARMED-STRANGE MESONS WITH THE RADIAL 

EXCITATION 



Before illustrating the strong decay of P-wave charmed-strange mesons with the radial excitation, we first introduce 
the category of the heavy flavor meson. 

In the heavy quark limit mg — > oo, the heavy quark plays a role of a static color source to interact with the light 
part within the heavy flavor hadron. Thus, the spin of the heavy quark sq can be separated from the total angular 
momentum J of the heavy flavor hadron. Furthermore, ji = Sg -I- L is a good quantum number, where Sq and L 
denote the spin of the light part of the heavy flavor hadron and the orbital angular momentum between the light part 
and the heavy quark, respectively. 

Thus, the heavy mesons can be grouped into doublets according to jf, which include je = ^ doublet (0^, 1^) with 

the orbital angular momentum L = 0, je = doublet (0+, 1+) and |^ doublet 2+) with L = 1. For L = 2 there 

exist (1~,2^) and (2^,3^) doublets with jf = | and | , respectively. As showing in Fig. [U the states existing in 
the doublets (0",1~), (0+,l"'') and (1"'",2+) are already filled with the observed charmed-strange mesons. Two 1+ 
states existing S and T are the mixture of two basis states l^Pi and I'^Pi [2^. l23j 




cos 6 sin 9 
— sin 6 cos 9 



where one takes the mixing ang le 6* = -tan-i V2 = -54.7° according to the estimate in the heavy quark limit. 

For P-wave charmed-strange mesons with the radial excitation discussed in this work, one also categorizes them as 
S = (0+, 1+) and T = (1+, 2+) doublets according to the above approach. Two 1+ states are the mixture of two 
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basis states 2^ Pi and 2^ Pi, which satisfy the below relation 




cos 6' sin 9' 
— sin 6' cos 6' 



2'Pi) 
2^ Pi) 



In this work, we approximately take 6' — 9 = —54.7°. 

For distinguishing P-wave states with and without the first radial excitation, one labels four P-wave states without 
the first radial excitation as 0"'"(S'), l^(S'), ^^{T) and 2+(T). Four P-wave states with the first radial excitation are 
named as 0+(S'*), 1+(S'*), l+(r*) and 2+(r*). 

If we set the upper limit of the masses of P-wave states with the first radial excitation as 3.04 GeV, the kinematically 
allowed decay modes of P-wave states with the first radial excitation are listed in Table [H In the following, the "^Pq 
model will be applied to calculate these strong decays in Table HI 



State 



1+(S*)/1+(T*) 



2+(T*) 



Decay modes 



0" -hO" 
1^ + 1" 
0+ + 1" 
l+(5)-f 0- 
l+(5') + r 

i+(r) -1-0" 

1+{T) + 1- 
2+ + 1" 



0" -f 1- 

I- + 1- 

0+ -hO" 

0+ + 1- 

l+[S) + Q- 
1+{T) + Q- 
l+(5)-f 1- 
1+{T) + V 

2+ + Q- 
2+ + 1- 



0" +Q- 

1- -hO- 
1- -hi- 
0+ + 1- 
i+(5) -hO- 
l+(S')^- r 
i+(r) + 0" 
i+{T) + r 
2+ + 0- 

2+ + 1- 



Decay channels 



D+K\ D°K+, Dtn^''> 
D'+K'\ D*''K'+ 

X 

L»i(2430)°7f+, Di{2A30)+K°, D,i(2460)+?7 



Di{2420)+K°, Di{2420)°K+ 

X 
X 



D+K*", D°K*+, D+<j) 
D*+K\ D*°K+, Dt+v 
D'+K*°, D*°K'+ 
DS(2400)+Jf°, DS(2400)''A'+, D:o(2317) + r? 

X 

L»i(2430)''A'+, i3i(2430)+A'", I>,i(2460)+77 
Di(2420)+K°, Di{2'i20)°K+ 

X 

X 

02(2460)+ K°, D|(2460)°is:+ 

X 



0+A'°, D°K+, 75+7,0 
D+K'^, D^K''+, Dt(t> 
D*+K°, D*°K+, D*+r) 
D'+K'\ D*''K*+ 

X 

Di(2430)''A:+, Di(2430)+A^ Dsi(2460) + ?7 

X 

I?i(2420)+A°, Di(2420)''A:+ 

X 

L»2(2460)+A'^ D2{2460)°K+ 



TABLE I: The relevant strong decay modes of P-wave charmed-strange mesons with the first radial excitation allowed by the 
conservation of the quantum number. Here "x" denotes that the decay modes are kinematically forbidden if setting the upper 
limit of the meisses of P-wave states with the first radial excitation as 3.04 GeV. Since the 1+ state in the (l+,2''") doublet 
decays into D*n via D-wave, it is very narrow and denoted as -Di(2420) The 1+ state in the (0+, 1+) doublet decays into 
D*Tv via S-wave. Hence, it is very broad and denoted as -Di(2430) 
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A. A review of the QPC model 

The '^Pq model, also known as the Quark Pair Creation (QPC) model, was firstly proposed by Micu in Ref. to 
calculate Okubo-Z weig -Iizuka (OZI) allowed strong decays of a meson. Later, this model was developed by the other 
theoretical groups |^ |^ 
decay of hadron (gI. |22|. |23|.12 



— I 

25, 



2C, 
2f , 



21 



and is successful when applied extensively to the calculation of the strong 

In the QPC model, the heavy flavor meson decay occurs through a quark-antiquark pair production from the 
vacuum, which is of the quantum number of the vacuum, i.e. 0++ For describing a strong decay process of 

the charmed-strange meson A{c{l)s{2)) —* B{c{l)q{3)) + C{s{2)q{4)), one writes out the S-matrix 

{BC\S\A) =1- i2TT5{Ef - E,){BC\T\A). (1) 

In the non-relativistic limit, the transition operator T is depicted as 

T = -37^(1to;1 - to|0 0) /dkg dk4(53(k3 + 

m •' 
Xy,„. (^^^) xlUn Vl' 4(k3) &i,(k4) , (2) 

where i and j denote the SU{3) color indices of the created quark and anti-quark, ipo^ ~ {uu + dd + ss)/V3 and 
cjq"* = -^<5q3Q4 {a = 1, 2, 3) correspond to flavor and color singlets, respectively. Xi^-m ^ triplet state of spin. 

y^mO^) = \M^Yem{(^k, (f'k) is the £th solid harmonic polynomial. 7 is a dimensionless constant which represents the 
strength of the quark pair creation from the vacuum and can be extracted by fltting the data. 
For the convenience of the calculation, one usually takes the mock state to depict the meson [13] 

\Ain'S+'LjM,,){K)) 

= V2E J2 {LMlSMs\JMj) 

Ml. Ms 

X J dkidk2<53 (K - ki - ka) *„LM, (ki,k2) 

xxFM.^''^''ki(ki)g-2(k2)), (3) 

which satisfies the normalization conditions {A(K)\A(K')) = 2E 6^(K — K'). Here, '^uLMl (ki,k2) is the spatial 
wave function describing the meson. 

Taking the center of the mass frame of the meson A: Ha = and Kg = — Kp = K, further we obtain a general 
expression of eq. H]) 

{BC\T\A) - ^SEaEbEc j J2 (1 1 - "^1 0){LaMl_^SaMs^\JaMjJ 
X {LbMl^SbMsAJbMjs){LcMlcScMsc\JcMjc) (<Pb <^c I^aVo^) 

X {X^SbMss xfcMsc \Xs\ms^ xWn)J^MLg',MLc ("^^ 

The color matrix element {u)]^ Lojj\Lo'^^ lu'^'^) = 1/3, which cancels out the factor 3 before 7 in cq. {f^s 'P'c^\f^A 'Po^) 
and (XsLa/s X^ScMs \x^SaMs Xi^) ^-re the flavor matrix element and the spin matrix element, respectively. Here, 

the spatial integral Im^^'^^ 



CZ%Lc = y"dkidk2dk3dk4 ^^(ki + k2)(53(k3 + k4)52(KB - ki - \ii^)5\Kc -^2- 

><*,*«L«M^^(kl,k3)*:^i^M^^(k2,k4)*„^L^A/^^(ki,k2)3^l™(^^2-^)^ (5) 

which describes the overlap of the initial meson [A) and the created pair with the two final mesons [B and C). 
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In this work, we use the simple harmonic oscihator (HO) wave function to represent the radial portions of the 
meson space wavefunction. The wave functions corresponding to the states with nL = IS", IP, 2P are respectively 



V'n=l,L=o(k) 
V'n=l,L = l(k) 
V'n=2,L = l(k) 



^3/2 



-i2 



exp 



3 7ri/4 

2^5/2 



3^r(k)cxp 



157ri/4 



(5-2k2i?2)3;r(k)cxp(-^), 



(6) 
(7) 
(8) 



which satisfy the normalization J ■(/;* ^(k)^„^i(k)(ik = 1. Here the solid harmonic polynomial 3^{"(k) = y ^ e_„j • k 

with e±i = (±l/-\/2, — «/v^, 0) and eq = (0,0,1). k = (m^kj — mjk.i) / (irii + rrij) is the relative momentum between 
the quark and the antiquark within a meson when considering the quark mass difference. These HO wave functions 
are relevant to the calculation of the strong decay of P-wave states with the first radial excitation. 
The helicity amplitude satisfies the relation 



{BC\T\A) = S^{Kb +Kc- Ka)M^^''-^^'-'b^''-'c . 
In terms of the partial wave amplitude, one obtains the partial decay width 



A/f2 



Af 



JL 



(9) 



(10) 



where |K| denotes the three momentum of the daughter mesons in the parent's center of mass frame. The partial 
wave amplitude A4'^^ is related to the helicity amplitude M^^-'a^^-'b^'-'c via the Jacob- Wick formula [35j 



BC) 



2Ja 



J J2 {LOJMjJJaMj, 



{JbMj, JcMj^ I JMjJM''''^ ''''B ''■'c (K) , 



(11) 



where J ~ Jb + Jc a-nd J a + Jp = Js + Jc + L. The calculation of transition amplitude using the '^Pq model involves 
two parameters: the strength of quark pair creation from vacuum 7 and R in the harmonic oscillator wavefunction. 
7 is an universal parameter which was already fixed from other channels as indicated in Ref. (sH]. The value of R 
is chosen to reproduce the Root Mean Square (RMS) radius obtained by solving the schrodinger equation with the 
linear potential. 

B. The partial wave amplitude of two-body strong decays of P-wave states with the first radial excitation 

With the preparation mentioned above, we obtain the partial wave amplitude of the strong decays of the P-wave 
states with the first radial excitation, which are listed in Tabic [D In Tabic |TT1 the concrete expression of the partial 

wave amplitude is given. The details of the spatial integral I^'^^'^j^ (K) are given in the appendix. 



State 



Decay channel 



Partial wave amplitude 



0+(S*) 



0" +0" 

i+(5) ^-o- 
l+(^) -f 0" 



M''° = ^VEAEBEc-y[2I± - Jo] 
= f^VEAEBEc-yilo - 2J±] 



=cos6'[- ^^EaEbEci{2II-^ - /, 
+ sin 



- sin 6* [ 
4- cos e [ 



fVEAEBEcj{llr, - 
-fVEAEBEc"/{2Poo' - 



00 I 



2a 
3 



VEAEBEc-yinS-l!^,)] 



Q- + r 



1- +0" 



1" + r 



0+ +0" 



l+{T) + 0- 



2+ +0" 



M 



cos e'[f 



VEAEBEc'y{2I±- lo)] 
+ sin [ - ^ VEaEbEcJ (2/± - Jo ) 
= cose'[^VEAEBEci{2I± + 2Io)] 
+ sme' EaEbEc-i{I± + lo)] 

7W 1" = cos 6*' [ I ,JI^EaEbEc-i (2/± - /o) ] 
+ sin e' [ 



2a 

3^3 ' 



-.^EaEbEc'){2I±- h) 

^^^EaEbEci{2I± + 2/o)] 
+ sme' [^^/EAEBEc'y{I± + Io)] 
X"' = cos [ - ^ ^VEAEBEc'y{2I± - Jo) 
M""^ = sine' [^VEAEBEcn{l± + lo)] 



cosO'lf 



\VEAEBEc-/{lSS+2f^^) 
+ sin e' [^^EaEbEc^{II-^ + JS)] 
= sin 6* cos 6»' [^^EaEbEci{2IIo^ + 21 
+ cos e sin [ - ^^EAEBEcy{llo - I- 
+ sin 6 sin 6*' [ - §VEaEbEc-/{ - J,?o - l': 
" = cos e cos [^^EaEbEcj (2 JoV ^ + 2JS 



r01\ 
'10 ) 



rOl 
'10 



- sin Osine'l^^ VEaEbEc-/ {llo - l': 
+ cos 6 sin 

cos e' 



- fVEAEBEc-yi 
[j^VEAEBEc^i^n 
- j^yEAEBEcijSll 



rOO 
J 00 



4/i'(5-6JoV 



'00 

+ 6Ji\,, 



rlO\ 



rOO 
'00 



+ 2J!, 



i+(r* 



0- + r 



1- +0" 



1" + 1" 

0+ +0" 

i+(S') + o- 
i+(r) + 0" 

2+ +0" 



= ^ sine' [% y EaEbEc-i{2I±- Iq)] 

+ cose'[~ ^VEAEBEc'y{2I±- lo)' 
^ -shie'[^y/EAEBEci{2I±+2h)] 

+ cos9'[^VEAEBEc-/il± + lo)] 
= - sin 6' [f yi ^/EZBi^7 (2J± - Jo) ] 

+ cose'[ - ^^EaEbEc'){2I± - Jo); 
= -sin6l'[^VS.4£s£c7(2J± +2Jo)] 
+ cos 6»' ^/EIEi£^7 ( J± + Jo) ] 
M"-'' = -sin^'f- ^^^EAEBEc'y{2I± - Jo) 
Al22 = cose' [2^^EZEi£^7(J± + Jo)] 



-sin^'tf 



^S) 



iVEAEBEc-YiC +21^ 
+ cos [|^/BZBi£^7(JoV' + nS)] 
= - sine sin 6*' [^^EaEbEcj{2I^o'' + 2llS)] 



- cos ecos9'[-^ VEaEbEci ( J^ - J?o' ) ] 

f VEaEbEc-J ( - Jc?g - JS + loo 



- sm t/ cos f 



+ cos 6 cos f 



\^VEaEbEcj{2I, 

- ^^EAEBEcy{l 1 

- fy/EAEBEc-fi^n 



00 



M^i = _ sin e' [^^^EAEBEcy{4lSS - 4J?,1 - 6J, 



1-1 

00 
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+ cos e' [ - ^^EAEBEcjiSlSS - - 71^,0 ' + 2/i?)] 




Q- +0- 


M"^ = ^^EAEBEc'y[I± + Io\ 




Q- + l- 


M'^ = j§^VEAEBEcj[l± + Jo] 




1- +0- 


M^^ = j§^VEAEBEc'y[l± + Jo] 


2+{T*) 




= ^^^EaEbEc^[2I± - /o] 




i+(s) + o- 


= cos e [-^^EaEbEci{4C + 6/S + 4/^7^ + 6/^)] 
+ sin e[^VEAEBEcj{3C + 7I?i + S/qV + 2/^)] 




1+(T) + 0^ 


= - sme[j^VEAEBEcj{^C + 6/?o' + 44"' + 6Jjg)] 
+ cose[^y EAEBEc-y{3C + 7/?o' + 3/oV' + 2/S)] 




2+ +0~ 


= ^,VEaEbEcj[C - /?o^ + 47' + 4/S] 



TABLE 11: The expression of the partial wave amplitude for the strong 
decays of P-wave states with the first radial excitation. Here a = 
2/VT8, —1/3 are for the strong decay involved in rj and rj' mesons re- 
spectively while a = l/y/3 is for the other strong decays, which are the 
result from the flavor matrix element. 



C. Numerical result 

The input parameters involved in the '^Pq model include the strength of quark pair creation from the vacuum, the 
R value in the HO wave function and the mass of the meson. One takes the strength of quark pair creation from 
the vacuum as 7 = 6.3 [36| . which is ^/96t: times larger than that adopted by the other theoretical groups 
The strength of ss creation satisfies 7^ = j/V^ If reproducing the reaHstic root mean square (RMS) radius by 
solving the schrodinger equation with the linear potential, one can obtain the value of R in the HO wave function 
[3^. The mass and the R value used in this work are shown in Table [1111 



State 


mass (MeV) [Ij 


R (GeV"') [36| 


D 


1869.62 (±) 1864.84(0) 


1.52 


Ds 


1968.49 (±) 


1.41 


D* 


2010.27 (±) 2006.97(0) 


1.85 


Dl 


2112.3 (±) 


1.69 


1)^(2400) 


2403 (±) 2352 (0) 


1.85 


D*o(2317) 


2317.8 (±) 


1.75 


-D.i(2460) 


2459.6 (±) 


1.92 


Di(2430) 


2427 (±) 2427 (0) 


2.00 


Di(2420) 


2423.4 (±) 2422.3 (0) 


2.00 


1)2(2460) 


2460.1 (±) 2461.1 (0) 


2.22 


K 


493.677(±) 497.614 (0) 


1.41 


K* 


891.66(±) 896.00 (0) 


2.08 


V 


547.853 


1.41 


n' 


957.66 


1.41 


</> 


1019.455 


2.08 



TABLE IIL The R value in the HO wave function and the mass relevant to the strong decays listed in Table H) Here (±) and 
(0) denote the charge of the meson. 



If the mass of the charmed-strange meson with 0"'"(5*) is 2.837 GeV predicted in Ref. [T3|, there only exists the 
decay channel 0^{S*) ^ 0~ -|- 0~, which is allowed by the phase space. In Fig. [21 we give the dependence of the 
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FIG. 2: The variation of the strong decay mode 0^(5*) — > 0^ + 0^ with Ra- 

partial decay widths of the strong decay of 0+(5*) state on the Ra- Here, Ra is the R value of the HO wave function 
of charmed-strange state with 0+(S'*). The minimum of the decay width around Ra = 1-7 GeV^^ in Fig. [2] is due to 
the node in the radial wave function of 0+(S'*). The total decay width of 0+(5*) charmed-strange meson is 108 McV 
with Ra = 2.8 GeV-^ 

In this work, we take the masses of and 1+(T*) charmed-strange mesons as 3.044 GeV, which is the 

experimental value of the mass of Z3s,/(3040)"'". Then, we calculate the decay of Dsj(3040)"'" under the two assumptions 
1"'"(5'*) and 1+(T*). In Figs. [3]and[4l we present the numerical results of the two charmed-strange mesons 
and l"'"(r'^). The dependence of the total decay width of on the Ra is shown in Fig. [6l Here. Ra denotes 

the R value in the HO wave function of £'sj(3040)"'". By comparing the calculated total decay width of Dsj(3040)+ 
with that of the Babar data, one finds that the total decay width (~ 204 MeV) of Dsj{304:0)^ obtained by the "^Pq 
model reaches up to the lower limit of the experimental width of I?sj(3040)"'" when taking Ra as 2.8 GeV~^. With 
increasing the Ra up to 3.5 GeV""'^, the total decay width is close to the central value of the experimental width 
of i3sj(3040)+. Thus, studying the decay width of Z?s,/(3040)+ under assignment shows that the first radial 

excitation of P-wave charmed-strange meson to Dsj(3040)+, i.e. is suitable. 

The result of the partial decay widths of l+(5'*) charmed-strange meson corresponding to Ra ~ 2.8 GeV~^ (see Fig. 
O indicates that 0" + 1" {D+K*°, D^K*+ and and 1" -I- 0" {D*+K°, D*°K+ and Dl+rf) are the dominant 

decay modes of I?sj (3040)"'", which further explain why Dsj(3040)"^ was firstly observed in D*K decay channel. An 
experimental search of I?sj(3040)+ in 0^ + 1^ channel {D'^K*^^ D^K*'^ and Dfcf)) is encouraged in terms of the 
ratio 

r(i+(5*)^o- + i-) 
r(i+(5*) ^ 1- + 0-) ~ ■ 

corresponding to Ra = 2.8 GeV~^. 




Under the assignment of 1^{T*) to _Dsj(3040), we can obtain the variation of the strong decays for l+(r*) 
0" + I- + 0-, 1- + 0+ + 0", l+(5) + 0", 1+(T) + 0", 2+ + 0" with the factor Ra {R value of the HO 
wavefunction of 1+(T*)) which is depicted in Fig. 31 Furthermore, the dependence of the total decay width on the 
Ra value is hsted in Fig. [T] The total decay width of D^j(3040)+ is about 33.8 MeV with Ra = 2.8 GeV~\ which 
is consistent with our knowledge, i.e., the 1^ state existing T doublet is of narrow width. In fact, the result of the 
decay of l"'"(r*) state further indicates that Z?sj (3040)+ can not be explained as l"'"(r*) charmed-strange meson. 

One also predicts that the partial decay widths corresponding to the decay channels 1+(T*) 0^ + 1^, 1^+0^, 1^ + 
1", 0+-H0-, l+(5)-hO-, 1+{T) + Q-, 2+ + are 9.8 X lO'^ MeV, 6.3 MeV, 13.0 McV, 10.1 McV, 9.9 x lO^^ MeV, 
3.5 MeV and 1.3 X 10~i McV, respectively These numerical results show that 1" -HO", 1" -fl", 0+ -f 0", 1+{T) + Q~ 
channels are important when searching l+(r*) state in experiment. 
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FIG. 3: The variation of the strong decays for (a) 1+(S'*) -> 0" + 1~, (b) 1+(S'*) -> 1~ + 0~, (c) -> 1" + 1", (d) 

1+(S'*) ^ 0+ +0", (e) ^ 1+(S) +0", (f) 1+(S'*) l+(r) + 0" and (g) 1+(S'*) ^ 2+ + 0" with the factor Ra of the 

HO wavefunction of 1"'"(S'*). Here the total partial decay width is labeled by "total" in diagrams. 



The dependence ofthe Strong decays 2+ (T*) ^ O' + O", Q-H-l", I'+O", I' + l", l+(S') + 0-, l+(T) + 0-, 2++0- 
on the factor Ra {R value of the HO wavefunction of 2+(T*)) is given in Fig. [5l Here, we take the mass of 2+(r*) 
as 3.157 GeV When taking Ra = 2.8 GeY-\ the total decay width of 2+{T*) is 87.9 MeV (see Fig. [T]), and the 
partial decay widths (see Fig. [5]) respectively corresponding to 2+(T*) ^ 0~ +0~, 0~ + 1~, 1~ +0~, 1~ + 1"'"(S') + 
0", l+(r) + 0", 2+ + 0" are 15.6 MeV, 0.49 MeV, 7.2 MeV, 49.3 MeV, 1.8 MeV, 13.2 MeV, 0.28 MeV, which show 
that 1^ + 1^, 0^ + 0^, 1^(T) + 0~ and 1^ + 0^ are key decay channels to find 2+(r*) charmed-strange meson. 



III. SUMMARY 



Stimulated by the newly observed charmed-strange meson Z?sj (3040)^^, we systemically study the two-body strong 
decays of P-wave charmed-strange mesons with the first radial excitation. 

Our numerical results show that Dsj(3040)+ can be categorized as 1+ state in 5 = (0+,l+) doublet well, i.e. 
-Dsj(3040)+ is the first radial excitation of Dsi(2460)"''. We suggest experimentalist to search _Dsj(3040)+ by 0^ + 1~ 
channel {D+K*^, D°K*+ and D+(j)). 

In the past six years, Babar and Belle experiments have made big progress in searching for charmed-strange mesons, 
which lets us believe that more charmed-strange mesons will be found in future experiment. If 1)57(3040)+ is the 
first radial excitation of Dsi(2460)"'", there must exist three partners of Dsj(3040)"'", which are the rest three P-wave 
charmed-strange mesons with the first radial excitation. In this work, we also study the strong decays of the rest three 
P-wave charmed-strange mesons with the first radial excitation. Our numerical result (see the presentation in the 
subsection of numerical result) will be helpful to instruct future experimental search of the remaining three P-wave 
charmed-strange mesons with the first radial excitation. 
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FIG. 4: The variation of the strong decays for (a) l+(r*) -> 0" + 1~, (b) l+(r*) 1" + 0", (c) l+(r*) ^ 1" + 1", (d) 
l+(r*) ^ 0+ +0", (e) l+(r*) ^ l+(S')+0", (f) l+(r*) ^ l+(r) + 0" and (g) l+(r*) ^ 2+ with the factor Ra of the 
HO wavefunction of l"'"(r*). Here the total partial decay width is labeled by "total" in diagrams. 



Note added. When this manuscript was completed, a work of Dsj(3040)"'" appeared (s^l- In this work, authors 
investigated the Ds mesons by a semi-classic flux tube model and explained Dsj (3040)+ as l+(j^ = 5 ). In our 
case, we calculated the strong decays of Ds,/(3040)+ with the assignment of the first radial excitation of Dsi(2460). 
By comparing the total decay width of Dsj(3040)+ obtained by the ^Pq model with the Babar data, we conclude 
that £'sj(3040)+ is the first radial excitation of Dsi(2460), which is consistent with the conclusion of the structure of 
Z?.j(3040)+ m. 
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Appendix 

According to the the spatial integral in eq. (O, one can categorize the strong decays of P-wave charmed-strange 
mesons with the first radial excitation into two groups: 2P IS" 15* and IP IP + IS. 

For the case of 2P IS + IS, the spatial integral ^^/^'^'^i/^ is simplified as Im'n' due to Ml^ = Ml^. = 0. 
According to the constraint from eq. (fTTj) . we take the direction of K along z axis: K ~ (0, 0, |K|). In the following, 
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